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Abstract 

We introduce (binary) Darboux transformation for general differential 
equation of the second order in two independent variables. We present a 
discrete version of the transformation for a 6-point difference scheme. The 
scheme is appropriate to solving a hyperbolic type initial-boundary value 
problem. We discuss several reductions and specifications of the transfor- 
mations as well as construction of other Darboux covariant schemes by 
means of existing ones. In particular we introduce a 10-point scheme which 
can be regarded as the discretization of self-adjoint hyperbolic equation. 

Integrable systems, Jonas transformations, Moutard transformations 



1 Introduction 

One can observe increasing role of difference equations over the past few decades. 
Primarily efforts were undertaken to discretize differential equations so that not 
to lose the properties (e.g. symmetries) that differential equations exhibit. It 
turned out quickly that difference equations in many aspects are richer and more 
fundamental than their continuous counterparts (many interesting structures 
disappeared under a continuum limit) and difference equations started to be 
something more than equations mimetic differential equations. In the present 
work we encounter the essential differences between discrete and continuous 
mathematical structures once more. 

The aim of this paper is to complete existing theory of Darboux transfor- 
mations (or better to say Darboux-Moutard transformations [1, 2, 3] and Jonas 
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transformations [4]) for differential equations and what more important to show 
the impact of the generahzation on the theory of Darboux-Moutard transfor- 
mations for difference equations. 

The main idea of this paper is to start systematic surveys that can free 
theory of integrable systems from their strong dependence of coordinate systems 
(parametrization of surfaces), desirable by many physicists author spoke to. It 
is especially important for the so called difference geometry [5] since due to 
results of the paper in the case of difference equations one can compensate lack 
of possibility of change of independent variables x = f{x, y), y = g{x, y). 

We recall that classical hmdamental transformation by Jonas [4, 6] regarded 
as the most general Darboux transformation acts on the conjugate nets in pro- 
jective space P" so in case when the net is two-dimensional it provide us with 
Darboux transformation for two-dimensional linear hyperbolic differential equa- 
tion in canonical form (from now on, unless otherwise stated, small letters de- 
note functions of real independent variables x and y and subscripts foregone by 
comma denote partial differentiation with respect to indicated variables) 

1p,^y+W1p,^+Zll),y+fll) = (1) 

and the transformation is nothing but the spatial part of (binary) Darboux- 

Backlund transformation for two-component KP hierarchy. When the net is 
more than two-dimensional the Jonas fundamental transformation provides us 
with Darboux transformation for the set of compatible equations of the form 
(1) and serves as spatial part of binary Darboux -Backlund transformation for 
multicomponent KP hierarchies (in other words yields Backlund transformation 
for n-wave interaction equations called sometimes Darboux equations see e.g. 
[7, 8]). 

The fundamental transformation has been successfully translated into the 
discrete language [9, 10, 11]. The discrete counterpart of conjugate nets are 
quadrilateral lattices in P" governed by system of equations of the type (un- 
less otherwise stated in the whole paper capital letters denote functions of two 
discrete variables m and n ((m, n) e Z^), and A„ denotes forward dif- 
ference operators ■= ^m+i.n — ^ and A„^' := 5" ,„.„+! — ^ while A_m 
and A_„ denotes backward difference operators A-m"^ ■= *m-i,n — * and 
A_„^ := ^^^,71-1 — note we identify = ^rn,n and in the whole paper we 
apply this convention) 

AmAn^ + AAm^ + BA„* + = (2) 

That is why in recent years notion of integrability of discrete (difference) equa- 
tions was often related to the planarity - 4-point schemes were the building 
blocks of the theory. 

Generalization to the so called quad-graphs (planar graphs), still objects 
with incorporated planarity, appeared only recently [12, 13, 14, 15]. 

In the present paper we show that planarity is not crucial from the point 
of view of integrable systems. It is remarkable but only example of more gen- 
eral theory with a 6-point difference scheme and 7-point self-adjoint difference 
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scheme as building blocks of the discrete theory of Darboux-Moutard transfor- 
mations. 

It turns out that the general differential equation of the second order in two 
independent variables 

{atjj^ +ctlj,y ) ,^ + (cV'.x +btp,y ) ,y +Wtlj,a: +ztp,y -ftjj = 

is covariant under a Darboux transformations (section 3) so conjugate nets are 
no longer of key importance. On the discrete level it reflects in the fact that 
one can generalize 4-point scheme to a 6-point scheme (see Fig. 1) 

A^m+2,n + B^m,n+2 + 2C*„+i,„+i + G*m+1,„ + i?*m,n+l = 

and quadrilateral lattices cease to be the master object of study in favour of 
triangular lattices with the 6-point scheme defined on them (section 4). 



m,n+2 



m,n+l m+l,n+l 



m,n m+l,n m+2,n 

• • • 



Figure 1: The 6-point scheme 
Moutard transformation for the 7-point self-adjoint scheme (see Fig. 2) 

^m+l,n-^m+l,n— 1 ~t~ ^m,n+l A^m-1,«+1 =.^iV, 

an example of equation given on a star (cross), has been derived in the paper 
[16] (see also subsection 7.2) and that is why we concentrate here mainly on 
the 6-point scheme. However, the present paper is thought to provide a brief 
overview on the topic of discretizations of 2D second order differential equations 
that are covariant under a Darboux transformations and the reader can find 
in the closing section 8 references to articles on integrable aspects of equations 
given on stars. At the moment we only underline that choice of difference scheme 
restricts sorts of initial-boundary value problems one can solve by means of the 
scheme. So it is important to indicate first on sort of initial-boundary conditions 
one would like to solve and then consider only the schemes that allow to solve 
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Figure 2 



7-point scheme 



the initial-boundary value problem. We start the paper with description of a 
well like initial-boundary value problem (section 2) we have in mind while 6- 
point scheme is considered. In turn the 7-point scheme is not proper to solve the 
initial-boundary value problem (as we know from numerical analysis schemes 
on stars are suitable to solve numerically Dirichlet boundary value problem for 
an elliptic differential equation, see e.g. [17]) but one can construct a proper 
scheme using the 7-point scheme (see subsection 7.3). 

In this paper we discuss also how the general Darboux transformation can 
be reduced or specified. We take the stand that introduction of novel termi- 
nology (such us specification) is necessary to discern procedures we deal with. 
We begin the discussion of reductions and specifications from the continuous 
case (section 5). Firstly, we consider the Moutard reduction (subsection 5.1) 
which is very classical construction [1] but to the best of our knowledge in full 
generality was given only recently [16]. Secondly, we discuss specifications (sub- 
section 5.2) and this part is (stands to reasons) new, specification to hitherto 
considered "conjugate" case or its elliptic counterpart are just examples of such 
procedure. Thirdly, we discuss two convenient gauge specifications (subsection 
5.3) the transformations can be written in. Finally we discuss reductions and 
specifications in the discrete case. 

We start from gauge specifications (subsection 6.1) and specifications (sub- 
section 6.2). We are not able to show reduction of the general 6-point scheme 
that leads to transformation of the Moutard type. Therefore we recall all results 
related to the topic first (section 7) and then all we are able to do is to introduce 
a 10-point scheme which is appropriate for solving defined in section 2 initial 
boundary value problem and can be regarded as a discretization of self-adjoint 
differential equation (section 7.3). 

We would like to stress once more that although we deal in the paper with 
linear equations only, the existence of Darboux- Moutard transformations makes 
this paper especially important for the theory of nonlinear integrable systems. 
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2 Well like initial-boundary value problem for 
6- point scheme 

In the present paper we pay special attention to difference schemes that allows 
to solve the following initial boundary value problem. We prescribe function 
\l'(TO,n) in the following points of the domain (see Fig. 3) 

• initial conditions 

{(m, n) € T\m + n = OV m + n = 1} 

• boundary conditions 

{(m, n) GT\{m = s — Ps An = Ps) \/{m = s—psAn = Ps), s = 2, 3,4, ...} 

where T denotes regular triangular lattice, Ps and Qs are functions N\{1} B s 
Ps such that Vs e N\{l}jJs < 
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Figure 3: Initial-boundary value problem. The initial values at points of two 
neighbouring straight line are given as well as two boundary conditions (black 
points). 

We concentrate in the paper mainly on the schemes that allow to find solution 
uniquely at least in the "upper half-plane" {(m, n) e T| m+n > 0} of the lattice. 
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For instance in the case of the 6-point scheme if the following conditions are 
satisfied 
1) 

V(m, n) e N X N, A„,„ 0, 7^ 0, 7^ 

2) 

Vs G N\{1} 

the matrices 

2C's_p^_2,pj. Bg-p^-2,ps • • • 

^S— Ps— 3,Ps + l -Pa-3,Ps + l -Bs-p,-3,ps + l 

■•• : 

: ■•• 

■•• As-q^-\.l^q^-l 2Cs~q^ + l^q^-l Bs-q, + l,q^-l 

■■■ ^s-gs,gs-2 '2Cs-q,,q,-2 

have non- vanishing determinant, then all the values at white points can be 
found uniquely. Similar result can be obtain for the 10-point scheme (49) with 
the only essential difference that the solution can be found uniquely only in 
upper half-plane. 

3 Darboux transformations for 2D second order 
differential equation 

It is a basic observation that map xj}* ^* given by 

■ip\x = Si)\a:+l3ip\y i)\y = -aip\^ --fip^y , (3) 

where a, j3, 7 and 5 are real functions (of independent variables x and y with 
an open, simply connected subset T) of as a domain) such that V(a;, y) G 
V, 7(5 — a/? ^ and + + (7 + (5)^ 7^ 0, is an invcrtiblc map between 
solution spaces of two differential equations of second order in two independent 
variables. Indeed, the compatibility condition of (3), which ensures existence of 
i/j* function, reads 

= 

(4) 

£* := adl + + (7 + 5)d:,dy + (a„ +5,y + (/?,„ +7,x )dy. 
Obviously i/i* satisfies equation of the same type but with bared coefficients 

As we shall see every second order equation in two independent variables 

C^xp = 

(6) 

:= adl + + 2cdxdy + (a,x +c,y +w)dx + {b,y +c,x +z)dy - f 
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can be transformed into the form (4) through a gauge transformation 

d ^ £,:= ^C^e (7) 

where ip and 6 are operators of multiplying by function (j){x,y) and 0{x,y) re- 
spectively, which we are going to determine now. Wc will call operator 
elementary transformable form of the second order differential operator. 

Indeed, request that operator C defined in eq. (7) is of the elementary 
transformable form £* (4) gives 

a = (l)ea I3 = (peh ^ + 6 = 24)60 (8) 

jC^e = (9) 



a,x+^,y = [{a9),x+{c9),y+a9,x+c6,y+w9](j), 
P,v+7,x = [{c9)„+{b9),y+c9„+b9,y+z9]cl). 



(10) 



On introducing auxiliary function p 



p:=^(^-7) (11) 

and treating eqs (8), (11) as the definition of fimctions a, (3, 7 and 5 and 
eliminating these functions from equations (10) one can rewrite equations (10) 
in the form 



(12) 



Function p exists provided that 

(t>d'9 - e{C^)^(t> = (13) 

Where {C^Y denotes the operator formally adjoint to the operator 

(£0^ := 9x{ad^ + cdy - w) + dy{bdy + cd^-z)-f (14) 

Taking into account (9) and (13) we come to the theorem 

Theorem 1 Gauge transformation (7) makes from arbitrary 2D second order 
operator U an operator in elementary transformable form iff 

£^61 = and (£^)V = (15) 

It turns out that just presented considerations lead to Darboux transformations 
for the 2D second order operator jC^ . Namely we have conclusion 

Conclusion 1 (Darboux transformations) We assume that 9 and (/> are 

class functions satisfying conditions (15), function p is given by formulae (12), 
r and s are arbitrary (of class C^) functions, and function d given by d := 
(p^ — (? + ab)(j) 9 obeys condition V(a;, y) € "D; d 7^ 0. Then the map ^ ■i^ given 
by 
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4>9 



p + c 

—a 



b 

p — c 




(16) 



is the map from solution space of equation ( 6) to the solution space of the equa- 
tion of the same form hut with the new coefficients 



= ° - - (17) 

£ := adl + hdl + 2cdxdy + {a„ +c,y +w)dx + {c,a; +b,y +z)dy - f ^ ' 

where the coefficients of (17) are related to coefficients of (6) by 

zz asr T hsr ^ csr 

"~<i' 

* = [t(f )>x +i(f ),. -(i)„ f ]r^ -z = \\{%y +i(f ),. +(§),. f p. (18) 

4 6— point scheme and its Darboux transforma- 
tions 

One can repeat considerations from the previous section in the discrete case. 
Indeed, starting from pair of equations 

A„§* = JA„** + /3A„«'* A„** = -aA„** - 7A„** (19) 

(whore the functions a, /3, 7 and 6 arc functions of discrete variables m and n) 
and writing down their compatibility condition 

am+l,n^'L+2,n + /^m,n+l*m,n+2 + (7m+l,ri + <5m,n+l)^'*n+l,ri+l~ 
{p-m-\-\,n + + 7rn+l,n 

(a + /3 + 7 + (5)** = 

(20) 

which is a 6-point scheme. One can ask if it is possible to transform the general 
6-point scheme of this type 

^*m+2,n + -B«'m,n+2 + 2C1'™+1,„+1 + G^^+x^n + i^^m.n+l = i^* (21) 

i^* = 

into the form (20) via a gauge transformation 

iF y^L:= ^L^e (22) 

only? The answer is positive and we will be calling equation of type (20) a 
6-point scheme in elementary transformable form. 
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Theorem 2 Gauge transformation (22) makes from the 6-point scheme (21) 
an operator in elementary transformable form L* iff the function O satisfies 
eg. (21) 

AQm+2,n + B<dm,n+2 + 2C0m+l,n+l + GQm+l.ji + H^m,n+l = (23) 

while function ^ is a solution of the equation formally adjoint to eq. (21) 



Am-2,n^m-2,n + B„in-2^n 



2Cm-l,n-l^m-l,n-l + 



(24) 



Then the functions a, j3, 7 and 5 in eq. (20) are given by 

^ — ^m — l,n^m— l,nQm+l,n 

P -^m,n— l^m,n — l®m,n+l 

7 ^ {Cm—l^n -Pm— l,n)^m— l,nQm,n+l 

5 = {Cm,n-l + -Fm,n-l)^m,n-10m+l,n 

where P is an auxiliary function defined by 



.(fa 



m+l,n+ 



iP) 



m.n — 



m — 1 , n m — 1 , n yjT, ^ n -(- 1 



cm 



m+l,n+l 



hi?$e™,„+i) 



(25) 



A_„($6to+i,„+iP) = ATO_i,„$TO-l,n0TO+l,n +^*0rn+2,n + 

Cm,n-l$m,n-10TO+l,n + C'^0m+l,n+l + Gi>0TO+l,rs 

We have the conclusion 

Conclusion 2 (Darboux transformations for the 6-point scheme) 

The map 1— »• # given by 



(26) 



A„(5*) 

{Cm,n—1 ~^ Pm.n — l)^m,n—1^7n-\-l,n ^'m,n—l^m.n — l^m.n-{-l 
-^m — l,n^m — l,n€)rn,+ l,n (-fm — l,n Cm—l,7i)^m — l.n^7n,n-\-l) 



A™ (f ) 

And) 



(27) 



where function Q satisfies eq. (21) while function $ is a solution of the eq. 
(24) and P is defined via equations (26), is the map from solution space of the 
equation (21) to the solution space of the equation of the same form with novel 
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(bared) coefficients related to the old ones via 



2C _ e„+2,„$^+i,„-i(c+p)„+i,„-i _|_ e^,„+2$m-i,,.+i(c-p)^-i,,.+i 



en. + 2,na'm + l,n-l(C+P)r^ + l,n-l+e,„ + 2,n^A 



I e^,„+i<i>^-i,„(c-p)^-i,„+e^+i,„$^-i,„A^_i,„ 

D ' 

H ^ e„.,n+2a-n.-l,n + l(C-P)„^i,„ + i + e^,„ + 2^B 

PS'„,„+i -D™,„ + i 

I e^+i,„$^,„-i (c+p)^,„-i+e„, „-iB^,n-i 

(28) 

where R and S are arbitrary non-vanishing functions, while D is a function 

given by 

D — [(-fm — 1 ,n C*m — l,n ) (^Pm,n—l''hCm,n—l )4"-Am — 1 ,n-Sm,n — l] ©m+l ,n ©m,n+l ^m — l,n^m,n — l 

and is assumed not to vanish on the whole lattice. 



5 Gauge specifications, specifications and reduc- 
tions 

One can reduce (or gauge or specify) the Darboux transformation in such a way 

that it maps between solution space of restricted class of equations. Assume 
first that matrix coefficients in eq. (16) and coefficients of its inverse obey the 
same linear constraint 

a'^^{p + c)+a'^H-a'^^a + a'^'^{p-c) = , . 

a^\p - c) - a^H + a'^^a + a^''{p + c) = ^'^^> 

where a'-' are given functions of x and y. Prom (29) we infer 

(a" + a22)p = 

and will discuss two cases p = and a^^ + a^^ = separately. 

5.1 Mouteird reduction (p = 0), reductions 

We have p = Q and (a^^ - a^^)c + a^'^b - a^^a = to be satisfied. The later 
constraint can be satisfied if one take a^^ = a^^, a^^ — and a^"^ = 0. To satisfy 
the equations (12) in the presence of condition p = it is enough to put (j) = 9 
andw = Q = z (i.e. demand that operator is formally self-adjoint). Moreover the 
functions r and s are no longer arbitrary, they must obey constraint ^ = const. 
As a result we obtain a transformation for formally self-adjoint equations which 
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arc usually refereed to (in the case a = = b, c = 2, s = r = ^9) as Moutard 
transformation. 

The procedure that impose the constraints on transformation's data (f) and 
9 {6 = (f) in the Moutard case) we call reduction of the transformation. 

5.2 Specifications a^^ + a^^ = 

The reduction is not the only procedure we have to our disposal. If we take 
a^^ = —a^^ we have the constraint 2a^^c + a^'^b — a^^a = to be satisfied. Two 

examples are 

a) a = = 6, a}^ = and c = 1 which is nothing but specification to 
"conjugate" case 

b) c = 0, a^^ = = a^^ (with option for deeper specification a = ±6) 

In those cases we only specify (specialize) the operator not affecting the 
transformation data. 



5.3 Gauge specifications, affine form, elementary trans- 
formable form 

The idea not to consider the operator itself but the equivalence classes with re- 
spect to the gauge goes back to Laplace and Darboux papers [18]. One can then 
develop theory in gauge independent language or choose a gauge appropriate 
to ones needs. We concentrate on the later procedure We take two arbitrary 
functions 9^, of class. Then operator 

:= (^01^90 

has coefficients 

(a^6^c^/») = ^V(a,&,c,L^^0) 



«;^=^T«^-(^0)^(^),.a-(^0)^(^)„c 

z<> = 9°cl,°z-{9or{p,yb-{9onp„c 
and operator adjoint to is 

(i9)t := 90Lf^o 

If in addition we define 

r = t, r = i e<> = i (30) 

then the form of the transformation remains unaltered 



{90cj>apO),^ = - (^9)2[c9(^),,+6.(^),^-^5^ 



e., (31) 
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(32) 



Wc indicate on two convenient gauges: 
1) Elementary transformable gauge 

L^e^ = {LS)U° = 

This gauge specification of the transformation to an elementary transformable 
form i.e., as we know from theorem 1, conditions 



r = 



+ Z0,y = 



hold. Moreover the operator {L^)^ is in elementary transformable form as well. 
The functions and are no longer arbitrary. To assure the transformed 
(bared) equation be in the elementary transformable form it is enough to put 
the functions and to be constant. 
2) Affine gauge is 



L^^o = 



so only 



r 







holds. Now it is enough to put to be a constant to obtain bared equation in 
affine form. Further convenient gauge specification is possible by putting 



Then the operator has the coefficients 

(a^, b^, c^, w^, z^, f^) = {6of{a, b, c, w, z, 0) 



6 Specifications, discrete case 
6.1 Gauge specifications 

In the discrete one can specify the gauge as well. Namely, we take two arbitrary 
functions G° $° Then operator 

has coefficients 

{A9,B9,C9,G9,H3,F3) = 

The two examples of convenient gauges are 

1) Affine gauge which we obtain demanding that 0° satisfies the equation 

L^e° = (33) 
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Then coefficients of operator obey constraint 

+ + 2CB + G^ + Hs -F^ = 

If one puts S = const then the constrain is preserved under the Darboux trans- 
formation. 

2) We demand that $° satisfies the equation 

(L^) t $° = (34) 
Then coefficients of operator obey constraint 

^m-l,n+l + ^m+l,n-l + ^C^ + + H^j^-^^ — „_|_]^ = 

If one puts R = const then the constrain is preserved under the Darboux trans- 
formation. 

If we apply the conditions (33) and (34) together and demand R = Q = S 
then we obtain transformation between elementary invertible forms of the 6- 
point scheme. 



6.2 Specifications, quadrilateral lattices, 3-point scheme 

A glance at transformation laws of the coefficients of the 6-point scheme (28) 
provide us with conclusions 

A) Both constraint A = 0, i? = and F = are preserved under the 
Darboux transformation (27). Constraint A — Q — B (or alternatively A = 
= F or B = = F) is specification to the celebrated 4-point scheme i.e. to 
quadrilateral lattice case subject of study of many papers (we confine ourselves 
to citing articles where Darboux transformations are considered) 

• Quadrilateral Lattices (Jonas fundamental transformations) [10, 11, 20, 

21, 22] 

• Circular Lattices and Quadratic Reductions (Ribaucour type transforma- 
tions) [25, 26, 27, 28] 

• Moutard type transformations [23, 24] 

• Symmetric (Goursat) type transformations [29, 30] 

Initial boundary value problem suitable for this scheme is no longer of the type 
mentioned in section 2. 

B) Constraint C = is not preserved under the Darboux transformation 
(27). So we have not got discretization of specification b) from the section 5.2. 

C) If we impose F = together with A = = B we obtain Darboux trans- 
formation for a 3-point scheme which corresponds to the continuous degenerated 
case -h 6^ -h = 0. 
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7 Discrete Moutard case 



The question arises can one find difference scheme which is appropriate to solve 
initial boundary conditions wc presented in section 2 and can be regarded as a 
discretization of Montard reduction? Wc were not able to find such a reduc- 
tion of 6-point scheme and we are not able to give a satisfactory non-existence 
theorem of such a 6-point scheme. Instead we constriict the Darboux covariant 
10-point scheme that can be regarded as a discrete Moutard equation and is 
suitable to solve initial boundary condition we described in section 2. 

In this section we firstly recall (for completeness of the paper) known results 
i.e. the discrete Moutard equation and its adjoint which were known only for the 
quadrilateral specification, so far. Secondly we recall self-adjoint 7-point scheme 
which is not proper to solve mentioned in the section 2 initial-boundary value 
problem. Finally we introduce a discretization of general Moutard reduction 
which is not just a reduction of 6-point scheme but arise from such reformulation 
of 7-point scheme that allows for solving the initial-boundary value problem. 

7.1 Discrete Mouteird equations 

On putting A = = B (quadrilateral specification), 2C = —F and G = H =: 
MF the equations (21), (23) and (24) take respectively form 

) (35) 

F(e„+i,„+i + 6) = G(e„+i,„ + e„,„+i) (36) 

^m,n-l (37) 

The crucial observation is: if the function © satisfies equation (36) then the 
function $ given by 

$ := -^(e™+i,„ + e„,„+i) (38) 
satisfies equation (37) [24] . If we put 

©m+l,n ©m,n-t-l 

$ 

then equations (26) will be automatically satisfied. If in addition we put S = Q 
then Darboux transformation (27) takes form 

A„e* = -ee„+i,„A™| ^""^^ 

which is the discrete Moutard transformation given by Nimmo and Schief [23]. 
Function ^ satisfies the equation 

^m+l,n+l + * = M{^m+l,n + ^m,n+l) (40) 
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where 



M = "" ',•1+^ (41) 

Equations (36) and (37) arc not the same so wc cannot say that wc have derived 
self-adjoint reduction. That is why we'd prefer to call the reduction Moutard 
reduction rather than self-adjoint reduction. 

7.2 7-point self-adjoint scheme and its 5-point specifica- 
tion 

To made the paper self-contained we derive in different manner the results 
contained in [16]. In the discrete case one can write 

A„iV* = 5A-mN* + p/^-nN* Ar^N* = -aA_mN' - 7A-nA^* (42) 

so instead of forward difference operators on the right side of equations (19) 
we have introduced backward difference operators in (42). The compatibility 
condition provide us with an elementary transformable 7-point scheme 

A„((JA_„iV* + /3A_„7V*) + A„(aA_„Af* + 7A_„iV*) = 

and function N satisfies an elementary transformable 7-point scheme as well. 
The question arises can one transform the general 7-point scheme 

to an elementary transformable 7-point scheme by a gauge transformation only? 
The answer is negative this time because apart from equations 

7 — Qm—l,n^7n—l,n^m,n—l- ^ '^ni,n — l^m,n—l^Tn—l,n^ 

a = -A<Pm-l,nO - Gm-l,n^m-l,n(^m,n-l (43) 

/? = —B$m,n-1& — y-m,n-l^m,n-l&m-l,n 

one has to satisfy three equations 

Q$&m-l,n + 1^^m,n-lOm-l,n = A$m-l,nO + C$rn-l,nOm,n-l 
H$0m,n-1 + C^m-l,n^m,n-l = ^^rra.n-lQ + I'^m,n-l6'm-l,n 

Happily enough the first two equations of can be satisfied by putting 

$ = e, g = A, n = B, c = v 

while the third one takes form 

*^m-t-l,n^m-t-l,n ~t~ -^^m— l,n ~1~ ^m,n-t-l^m,n-t-l ~t~ ^^m,ri— l~t~ (44) 
^m+l,n9m+l,n—l ~l~ ^m,n-t-l^m— l,n-t-l ~ ^9 

As a result we obtain Darboux transformation for 7-point self-adjoint scheme 
(c.f. [16]) 
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Theorem 3 Map N i—> N given by 







A„(iV0)_ 





m 
A ^ 



is a map from solution space of equation 

Cm+l,n-^m+l,n-l + Cm,n+l^m-'L,n+l = N, 

to solution space of equaMon 

Nm,n+1 + B Nm,n-1 + 
^m+l,n-^m+l,n— 1 ~t~ ^m,7i+l-^m— l,7i+l ~ !F N ^ 

where Q is a solution of equation (46) and the new (bared) fields are related to 
the old ones by 



(45) 
(46) 

(47) 



c 

F 



+l,n 



em-l.n-lPm-l, 

0{A 



-iV ' 



S. 



m,n+l 



(48) 



' ' ~Tn+l,n ~Tn — l,n 



'm,n+l 



-+ 



n 1 \p 1 



and where V := 6'^;B + 0m-i,nCA + Om,n-iCB . 

Clearly the scheme above admits specification C = (alternatively one can 
put A — QorB = Q) and as result we obtain specification to 5-point self-adjoint 
scheme (c.f. [16]). 

7.3 Discrete Mouteird type transformation, a 10-point scheme 

The idea is that in the Darboux transformation for the 7-point scheme we con- 
sidered in the previous subsection one can replace field A'^ with ^ m+i.n + '^ m.n+i 
without affecting leading terms of continuum limit. As a result we get the 10- 
point scheme 



[A + S)* + {Am+l,n + B„,„+i)*„+i,„+i-|- 

A}^m-l,n+l + B'^m+l^n-l + ^m+2,n + Bm,n+i^ m,n+2 + 

Cm+l,n^ m+2,n-l + Cm,n+l'^ m-l,n+2 + 

(Crre+l.n — F)^m+l,n + {Cm,n+1 — -P')*m,n+1 = 



(49) 



which is no longer self-adjoint but if the function O satisfies the equation 



(49) 



+ l)€)m+l,n+l + 
1 + B<dm+l,n-l + Am+l,n'dm+2,n + Bm,n+l^m,n+2- 
Cm+l,n^m+2,n-l + ^^,71+10771-1,71+2 + 

+ (^771+1,71 — -P')0r71+l,7l + (^771,71+1 " -P') ©771,71+1 = 



(50) 
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Figure 4: 10-point scheme 



then the function 

^ = 0m+l,n + Qm,n+1 

satisfies the equation formally adjoint to eq. (49) 

{A + B)<^ + {A l^m+l,n-l + 

^m — l,n+ 

Cm-l,n+l^rn-2,n+l + Crn+l,n-l'^m+l,n-2 + 

+ {C — Fm-l,n)^m-l,n + (C — -Fm,ra-l)^m,n-l = 



(51) 



The scheme is proper for solving the initial-boundary value problem we discussed 
in section 2. 

We finally receive transformation 



i)(H„, . 1. „-(-)„ 



i)) 



(e+e„ 



[A{em+l,n + Qm,n+l) + 0(0 + e„-^i,„-i)] A 



1'rri + l,,t+1'm.,t + l 

e™+i.„+e„,.„+i 

c(e + e„+i,„_i)A_„|^±i^±|^ 



(52) 



>+e^,„+i)) 



(0-Ke„+i,„_i) - 

— [B{Qm+l,n + ^m,n+l) + C(6 + 9m-l,n+l)] A_„^ 



c(e + e„_i,„+i)A_„|^±i^±|^ 
8 Conclusions and perspectives 

We end the paper with some conclusions and propositions for further develop- 
ment. 

The generalization of results of this paper to multi-dimension can bo given 
where set of difference operators A^, i — 1, N is replaced by a set operators 
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Di, i = 1, N that are linear commuting combinations of A^. The same reason- 
ing holds as far as Laplace transformations are concerned. The 6-point scheme 
admit factorization [{M-i_Ti+NiT2+Xi){M2Ti+N2T2+X2)+H]il; = (where Ti 
and T2 arc shifts in m and n direction respectively, so it can be used to develop 
theory of Laplace transformations for difference equations [31, 32, 33, 34, 24, 39]. 

We would like to mention that on this level of generality showing (/-difference 
analogue discrete schemes wc just introduced is straightforward [35]. 

Four propositions of further developments are in order. First, generalization 
of quadratic and symmetric reductions to 6-point scheme (it will be given in 
forthcoming papers). Second, simple idea by professor Decio Levi that not 
only operators A„ A-^ A_„ can be of importance e.g. one can try to 
use Ti — T2 and T1T2 — 1 operators instead. Third, to investigate the role of 
just presented transformations in the difference geometry. Finally and the most 
importantly, deriving hierarchies of nonlinear equations associated with all the 
equations presented in the paper. 

Since wc concentrated here on discretizations of hyperbolic equations we 
end the article with three comments on star schemes - proper discretizations 
of elliptic differential equations. Firstly, let us notice that if we substitute 
ip = (j) + (f'm.+i.n + 0m,n+i i^to thc 6-point schcmc wc will obtain star-like 
difference scheme, which can serve as star like discretization of eq. (6) and 
which needs further studies. 

Secondly, it is remarkable that star like (or cross like) operators such as 
7-point scheme or 5-point scheme appeared in the integrable literature occa- 
sionally [36, 32, 37, 33, 38, 39]. But almost none of the results were used to 
obtain solutions of nonlinear integrable systems. The only exceptions are works 
concerning discrete time Toda chains [40, 41, 42, 43, 44, 45] which are a non- 
linear 5-point scheme itself and the work [16] the result of which were used to 
obtain solutions of generalization of Toda chain to two dimensional lattice [46] . 
In other words theory of integrable difference systems based on the schemes 
other than 4-point schemes is still in its infancy. 

Finally, the relationship of integrable systems on quad-graphs with equations 
on stars the discrete time Toda type lattices are established [12, 13, 14, 15] 
and we would like to refer to the relationship as to sub-lattice approach. It is 
not clear under what circumstances the sub-lattice approach does not destroy 
integrability. However in the paper [47] it was proved that for discrete lattices 
governed by Moutard equation integrability features like existence of Darboux 
transformations are inherited by a sub-lattice. 
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